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H ' Abstract 

The temperature dependence of the liquid-drop fission barrier is considered, the 
critical temperature for the liquid-gas phase transition in nuclear matter being a 
parameter. Experimental and calculated data on the fission probability are com- 
pared for highly excited 188 Os. The calculations have been made in the framework 
of the statistical model. It is concluded that the critical temperature for the nuclear 
liquid-gas phase transition is higher than 16 MeV. 
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O i 1 INTRODUCTION 

The critical temperature for the liquid-gas phase transition is a crucial charac- 
teristic related to the nuclear equation of state. There are many calculations 
of I , for finite nuclei. In [I f2f5]H|5] . it is done by using a Skyrme effective 
interaction and the thermal Hartree-Fock theory. The values of T c were found 
to be in the range 10-20 MeV depending upon the chosen interaction param- 
eters and the details of the model. In Ref. [5117] the thermostatic properties of 
nuclei are considered employing the semi-classical nuclear model, based on the 
Seyler-Blanchard interaction. The value of critical temperature is estimated 
to be T c =16,66 MeV. 

As the temperature of a nucleus increases, the surface tension decreases and 
then vanishes at T c . For temperatures below critical, two distinct nuclear 
phases coexist - liquid and gas. Beyond T c there is not two phase equilibrium, 
only nuclear vapor exists. 

The main source of the experimental information for T c is the yield of interme- 
diate mass fragments. In some statistical models of nuclear multi-fragmentation 
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the shape of the IMF charge distribution, Y(Z), is sensitive to the ratio T/T c . 
It was noted in the earlier papers that the fragment charge distribution is well 
described by the power law, Y(Z) ~ Z~ T [8], as predicted by the classical 
Fisher droplet model for the vicinity of the critical point [Sj. In [8J the criti- 
cal temperature was estimated to be ~ 5 MeV simply from the fact that the 
IMF mass distribution is well described by a power law for the collision of p 
(80-350 GeV) with Kr and Xe. In the paper [10] the experimental data were 
gathered for different colliding systems to get the temperature dependence of 
the power law exponent. The temperature was derived from the inverse slope 
of the fragment energy spectra in the range of the high-energy tail. The min- 
imal value of r was obtained at T = 11-12 MeV, which was claimed as T c . 
The later data smeared out this minimum. Moreover, it became clear that the 
"slope" temperature for fragments does not coincide with the thermodynamic 
one, which is significantly smaller. A more sophisticated use of Fisher's model 
has been made in [TT]. The model is modified by including the Coulomb en- 
ergy release, when a particle moves from the liquid to the vapor. The data 
for multi-fragmentation in 7r (8 GeV/c) + Au collisions were analyzed. The 
extracted critical temperature was (6.7±0.2) MeV. The same analysis tech- 
nique was applied for collisions of Au, La, Kr (at 1.0 GeV per nucleon) with 
a carbon target [12J. The extracted values of T c are (7.6±0.2), (7.8±0.2) and 
(8.1±0.2) MeV respectively. 

Significantly higher critical temperature, 16.6±0.86 MeV, was obtained in [13] 
by semi-empirical analysis of the data for the "limiting temperatures" of frag- 
menting systems. Authors of Ref. [13] interpreted the obtained value as T c for 
the symmetric nuclear matter. 

Having in mind the shortcomings of Fisher's model [H|TT5] we have estimated 
the nuclear critical temperature in the framework of the statistical multi- 
fragmentation model, SMM [IB]. This model describes well the different prop- 
erties of thermal disintegration of target spectators produced in collisions of 
relativistic light ions. The intermediate mass fragment (IMF) yield depends 
on the contribution of the surface free energy to the entropy of a given fi- 
nal state. The surface tension coefficient of hot nuclei depends on the critical 
temperature. The comparison of the measured and calculated IMF charge 
yields provides a way to estimate T c . It was found from the analysis of the 
fragment charge distributions for the p(8.1GeV)+Au reaction that T c = (20 
± 3) MeV [T7]. In the next paper by the FASA collaboration [18] the value 
T c = (17 ± 2) MeV was obtained from an analysis the same data using a 
slightly different separation of the events. 

Thus, the different experimental estimations of the critical temperature from 
fragmentation data are very controversial. This is a reason to look for other 
observables which are sensitive to the critical temperature for the liquid-gas 
phase transition. It was suggested in Ref. [TS] to analyze the temperature 
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dependence of the fission probability to estimate T c . Note, that Silva et al. 
[20] explains why the power law used in the Fisher droplet model gives a 
spurious value for T c . 



2 TEMPERATURE DEPENDENCE OF FISSION BARRIER 



The fissility of heavy nuclei is determined by the ratio of the Coulomb and 
surface free energies: the larger the ratio, the smaller the fission barrier. As the 
temperature approaches the critical one from below, the surface tension (and 
surface energy) gradually decreases, and the fission barrier becomes lower. 
Thus, the measurement of fission probabilities for different excitation energies 
allows an estimate of how far the system is from the critical point. Temperature 
effects in the fission barrier have been considered in a number of theoretical 
studies based on different models (see e.g. [Tl[2ll[2 2 23 24 25|26] . The effect is 
so large for hot nuclei that the barrier vanishes, in fact, at temperatures of 
4-6 MeV for critical temperature T c in the range 15-18 MeV. 

In terms of the standard liquid-drop conventions [27J, the fission barrier can 
be represented as a function of temperature by the following relation: 

B f (T, T s ) = E S (T S ) - E° S (T) + E C (T S ) - E° C (T) = 

E° s (T)[(B s -l) + 2x(T).(B c -l)] (1) 



Here B s is the surface (free) energy at the saddle point in units of surface 
energy .E°(T) of a spherical drop; B c is the Coulomb energy at the saddle 
deformation in units of Coulomb energy E°(T) of the spherical nucleus; T s 
and Tare temperatures for the saddle and ground state configurations. For 
the surface energy and the fissility parameter x(T), one can write [2T] : 



E° S (T) = E° s (0) 



a(0) 



P(0) 
P{T) 



1 2 /3 



x(T) 



E° C (T) 



x(0) 



p(0)a(T) 



(2) 



where cx(T) and p(T) are the surface tension and the mean nuclear density for 
a given temperature. Equation (JT]) can be written as: 



B f (T,T s ) = B f (T s ) + AB f (3) 

where AB f = E®(T S ) — E®(T) + E®(T S ) — E®(T). Here B f (T s ) is fission barrier 
calculated under assumption that T s = T. In that case the values B s and B c 
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are determined by the deformation at the saddle point, which depends on the 
fissility parameter x(T). These quantities were tabulated by Nix [27] for the 
full range of the fissility parameter. The value of ABf is determined by the 
surface and Coulomb energies of a spherical drop, and can be easily calculated. 
For er(T) we use the approximation: 



a{T) = a(0) 



T 2 — T 2 

c 

T 2 + T 2 



5/4 



(4) 



This equation was obtained in Ref. [28] devoted to the consideration of ther- 
modynamic properties of a plane interface between liquid and gaseous phases 
of nuclear matter in equilibrium. This parameterization is successfully used 
by the SMM for describing the multi-fragment decay of hot nuclei. Figure 1 
shows the different approximations used in the literature for the surface ten- 
sion coefficient as a function of T/T c . 




T/Tc 



Figure 1. The calculated coefficient of the surface tension as a function of T/T c : 
lines 1 and 2 are obtained according to eq. (TJJ) and [5]), lines 3, 4 are for linear and 
quadratic parameterizations of c(T). The symbols are taken from Ref.[l]. 

Curve number 2 was calculated in framework of semi-classical model, based 
on the Seyler-Blanchard interaction [7J. An analytical expression for cx(T) 
obtained in this paper is the following: 



a{T) = a(0) 



T 

1 + L5 s 



T, 



1.5 



(5) 



Two other parameterization of er(T) are also presented: linear ~ (1-T/T C ), 
which is used in the analysis with the Fisher droplet model pdlfT2] . and 
quadratic ~ (1-T/T C ) 2 [29]. 
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In accordance with [27], the expressions for Eg(0) and x(0) are taken to be 



Z 2 1 A 

E°(0) = 17, 9439 7 • A 2/3 MeV, x(0) = '—. 

sK ' ' ' ' v ' 5O.887' 



(6) 



where 7 = 1 — 1.7826 — ^— . Sauer et al. pQ investigated the thermal proper- 
ties of nuclei by using the Hartree-Fock approximation with the Skyrme force. 
The equation of state was obtained, which gives the critical temperature T c 
18 MeV for finite nuclei. The temperature dependence of the mean nuclear 
density was found to be p{T) =p(0)(l - «T 2 ), where a = 1.26 ■ 1CT 3 MeV~ 2 . 
In the following we shall use this finding for p(T). 
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Figure 2. Relative value of fissility parameter, calculated for Z88 0s as a function of 
relative temperature for different parameterization of surface tension. Meaning of 
the lines is explained in caption of Fig.l. 

Figure 2 shows the relative values of the fissility parameter x{T) for ' 88 0s 
calculated as a function of reduced temperatureT/T c . This nucleus has been 
chosen since the results can be compared with well known experimental data 
for this nucleus [30J . The calculations are performed for the different versions 
of cr(T) mentioned above. A drastic change of nuclear fissility is expected even 
halfway to the critical point. 



Figure 3 displays the calculated value of the liquid-drop fission barrier for 
188 0s as a function of relative temperature. Virtually, the barrier vanishes 
for T > 0.4T C if the surface tension is taken according to (j4j) and (joT) . For 
the linear and quadratic approximations of a (T) the reduction of the fission 
barrier with temperature is much faster. 
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Figure 3. Temperature dependence of liquid - drop fission barrier for Os. The 
meaning of the lines is explained in Fig.l caption. 

3 THE ESTIMATION OF FISSION PROBABILITY 

In this chapter we analyze the experimental data on the fission probability 
of 188 Os, produced in collisions 4 ife+ 184 W |30j . The excitation energy of the 
compound nucleus created at the highest beam energy is 117 MeV. The shell 
and pairing effects are predicted to disappear for such a hot nucleus; there- 
fore the fission barrier is expected to be the liquid-drop one. This barrier is 
temperature dependent. Comparison of the measured and model calculated 
fission probabilities provides a way to estimate the critical temperature T c . 

Experimentally the fission probability Wf can be found from the measured 
fission cross section of. 



where <7r is total reaction cross-section. The main decay mode of the com- 
pound nucleus in 4 i?e+ 184 W collisions is the sequential emission of neutrons. 
For the highest excitation energy the mean number of emitted neutrons is 
11-12. The mean fission probability during a neutron cascade oix steps can be 
calculated by the following equation: 



The ratio Tf/T tot is the relative fission width for the i-step of the cascade. 
According to the statistical model [31] the value of T/is calculated as 



W f = a f /a R , 



(7) 
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Ui—Bfi 

Here U is the thermal part of excitation energy E*, p{U) is the level density, 
the index s is used for the saddle configuration. It is natural to use in (j3J) the 
temperature dependent fission barrier as has been done in [21.22 .23, 24.2 5|26] . 
The problem was considered also in [32]. The neutron width is given by the 
following equation [33] : 



r n (E*,I l ) = 2 ^ + ^™ n f a^EjpiiUi-Brt-EjEndEn (10) 



Here B ni , E n , S n are binding, kinetic energies and spin of the neutron, a n (E n ) 
is the neutron capture cross-section for the inverse reaction. The contribution 
of charged particle evaporation is on the level of several percent of T to tai- 
Nevertheless it has been taken into account. For level density p(U) the Fermi- 
gas model is used. 

Figure 4 presents the comparison of the data for fissility of 188 Os as a func- 
tion of excitation energy [20] with calculations under the assumption that the 
surface tension is described by eq. (4). The critical temperature is a parame- 
ter that can be found from the best fit. It is done for the highest excitation 
energy available, where the temperature dependence of the fission barrier is 
more prominent. 

The result is demonstrated in Fig. 5. Different calculations are presented, which 
have been done using all the parameterization of the surface tension mentioned 
above. It seems clear that the linear and quadratic approximations for cx(T) 
should be excluded as unrealistic. Fission probabilities, calculated with egs. 
(4) and (5) fall rather fast with increasing the critical temperature. They are 
crossing the experimental band giving the following values of critical temper- 
ature: T c w (23.5 ± 2.5) MeV in the first case, and at T c w (17.5 ± 1.5) 
MeV for the use of eq.(5). This is in accordance with the value of the critical 
temperature obtained by the FASA collaboration from multi-fragmentation 
data. These values are only slightly changed when the shell effect is taken into 
account for the last steps of the neutron cascade. 
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Figure 4. Fission probability of Os as a function of the excitation energy: dots 
are data [30J, curves are calculated assuming different values of critical temperature. 
Surface tension is taken according to (4). 
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Figure 5. Fission probabilities for Os at excitation energy 117 MeV. The cal- 
culated values (lines) are given as a function of the assumed critical temperature. 
Different parameterizations of surface tension are used (see Fig.l). The experimental 
value is shown by the horizontal band. 

4 CONCLUSION 

Critical temperature for the nuclear liquid-gas phase transition has been esti- 
mated from the fission probability of the highly excited nucleus 188 Os. Analysis 
is made under different assumptions about the temperature dependence of nu- 
clear surface tension. The results presented here provide strong support for 
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the value T c > 16 MeV. 
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